Implications of the recent high statistics determination of the pion electromagnetic 

form factor in the timelike region 
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The recently evaluated two-pion contribution to the muon g — 2 and the phase of the pion elec- 
tromagnetic form factor in the elastic region, known from tttt scattering by Fermi- Watson theorem, 
are exploited by analytic techniques for finding correlations between the coefficients of the Taylor 
expansion at t = and the values of the form factor at several points in the spacelike region. We 
do not use specific parametrizations and the results are fully independent of the unknown phase 
in the inelastic region. Using for instance, from recent determinations, (r^) = (0.435 ± 0.005) fm'^ 
and F(-1.6GeV2) = Q.243,tllfi, we obtain the allowed ranges 3.75 GeV^* < c < 3.98 GeV^"' and 
9.91 GeV~® ^ d < 10.46 GeV~® for the curvature and the next Taylor coefficient, with a strong 
correlation between them. We also predict a large region in the complex plane where the form 
factor cannot have zeros. 

PACS numbers: 11.55.Fv, 13.40.Gp, 25.80.Dj 



I. INTRODUCTION 

The pion electromagnetic form factor F{t), defined by 
the matrix element 



{7r+ip%r'^\n+ip)) = ip + p'),Fit) 



(1) 



where q — p — p' and t = q'^, plays a central role in 
strong interaction dynamics. From the general principles 
of quantum field theory, it follows that F{t) is normal- 
ized to F{0) — 1 and is a real analytic function in the 
i-plane cut along the real axis from the unitarity thresh- 
old t+ = \Ml to infinity. At low energies its proper- 
ties are described by chiral perturbation theory (ChPT), 
the low energy effective theory of the strong interactions 
P, Q , calculations of the pion form factor being available 
in ChPT up to two loops Lattice gauge theory 

has recently become another useful tool for the calcula- 
tion of the form factor at low energies Q. On the other 
hand, perturbative QCD predicts the behavior at large 
momenta along the spacelike axis, where = — i > 0. 
The leading order (LO) asymptotic term is [l|-[T2l| 



167r j^>,(Q^) 
Q2 



(2) 



where F.,^ is the pion decay constant and oigiCf'') — 
47r/[91n((5^/A^)] is the running strong coupling to one 
loop. Next-to-leading-order corrections to ^ were cal- 
culated by various groups As discussed, for in- 
stance in lU [3 ; transition to the perturbative 
QCD regime seems to occur quite slowly in this case. 

The experimental information available on the pion 
form factor is very rich. This quantity was measured 
at spacelike values > with increasing precision 



from electron-pion scattering and pion electroproduction 
from nucleons [l^-fl^. On the timelike cut, where the 
form factor is complex, the Fermi- Watson theorem im- 
plies that in the elastic region its phase is equal to the 
phase-shift of the P-wave of the tttt amplitude, calculated 
recently with precision using Roy equations and fixed-t 
dispersion relations (23 - 1261 . The modulus has been mea- 



sured from the cross section of e 



by sev- 



eral groups in the past [27|- [36l[, a nd more recently to 
high accuracy by the BABAR [33 and KLOE [H, IH 
collaborations. These data have been used for an accu- 
rate evaluation of the two-pion contribution to the muon 
anomalous magnetic moment [40l . IZlj . 

The constraints imposed on the pion form factor by 
analyticity and unitarity have been exploited in many 
works (the list [42|-[63| covers only partly a very rich 
literature). Different analytic representations , either as 
standard dispersion relations [T3|j phase (Omnes-type) 
[A H 113, m US mi or modulus ^ representa- 
tions, as well as expansions based on conformal map- 
pings [3 m, [Ml or Fade- type approximants 64] , have 
been constructed in order to correlate the low- and high- 
energy properties of the form factor. Of special interest 
is the issue of the zeros of the form factor, investigated by 
means of dispersive sum-rules [3 ES-SS lE^I or by the 
more powerful techni ques of analytic optimization the- 
ory [42, 113, . In (6ll - l63l . l66j | similar functional-analytic 
techniques were applied for deriving bounds on the ex- 
pansion coefficients at i = 0, from an weighted integral of 
the modulus squared along the cut, known from unitarity 
and dispersion relations for a related QCD correlator. 

In the present paper we address the same problem, i.e. 
to find constraints on the coefficients appearing in the 



2 



Taylor expansion 



F{t) = l + -{rl)t + ct^ + dt^ 
6 



(3) 



from a well-defined input on the timelike axis, and also 
include information coming from high precision experi- 
ments that measure the form factor in the spacelike re- 
gion. We also consider the problem of the zeros, and 
obtain a region in the complex i-plane where zeros are 
excluded. The main reason of revisiting the problem is 
the recent high statistics measurement of the modulus 
£(t)Jon the unitarity cut by BABAR [37] and KLOE 
38l |39j experiments. As we will show, this information 
leads to stringent constraints, of a remarkable level for a 
prediction independent of any specific parametrization. 

We apply a technique discussed in [6ll, [6^, which 
makes use of information on both the phase and mod- 
ulus, and was shown recently [H, [zO] to place stringent 
bounds on the Kit weak form factors. As first input we 
use the Fermi- Watson theorem, according to which one 
has, modulo tt, 



Aig[F{t + ie)] = 5l{t), t+<t<h 



(4) 



where 5\{t) is the phase-shift of the P-wave of tttt elastic 
scattering and tin the first inelastic threshold. As dis- 
cussed previously [H, , inelasticity in the case of the 
pion vector form factor is negligible below the opening of 
ITU) channel, so we take = {M.,^ + M^)"^. Below this 
energy, the phase Sl(t) is known with precision from Roy 
equations and fixed-t dispersion relations for tttt scatter- 
ing [Hi^. 

We also include information on the modulus, generi- 
cally expressed by an integral relation 



1 



dtp{t)\Fit)\^ < /, 



(5) 



where p{t) is a positive definite weight in the region of 
integration and / is a known quantity. Actually, ([5]) does 
not fully account for the present information on 
indeed, except for a small region near the threshold t-|- — 
4M^, the modulus is measured also below the inelastic 
threshold tin, i.e \F{t)\ is measured more or less point- 
wise, at every t, not only in averaged form as in In 
principle, the accurate knowledge of the phase and mod- 
ulus on a region on the unitarity cut is sufficient to pin 
down the form factor everywhere due to analyticity. In 
practice, however, due to the well-known "instability" of 
analytic continuation, the uncertainties, however small, 
lead to solutions which are very different at points outside 
the original data interval. Therefore, we do not proceed 
by constructing parametrizations of the form factor on 
the timelike axis, but consider instead the global class 
of functions compatible with the adopted input, and de- 
rive constraints on various quantities of interest from this 
class of functions. As we shall see, even the input ^ 
leads to quite strong constraints on the properties of the 
form factor near t = and in the complex plane. Thus 



the chosen method is fully justified by the results that 
have been obtained. 

A further open point is the choice of the weight p(t) in 
([5]). In principle, a large class of positive weights, lead- 
ing to a convergent integral for \F{t) \ compatible with the 
asymptotic behavior ([2]), can be adopted. The optimal 
procedure is to vary p(t) over a suitable admissible class 
and take the best result. This approach will be investi- 
gated in a future work. In the present paper we make 
the particular choice that corresponds to the two-pion 
contribution to the muon g — 2, when the weight p{t) has 
the form 



K{t) = f du- 
Jo 1 



t7/2 

(1 - u)u^ 



K{t), 



l-u + M2u2/t' 



(6) 



and the right-hand side (rhs) of JS]) is the two-pion con- 
tribution to the muon anomaly in the range t > tin, 



(7) 



The practical motivation of this particular choice is that 
an accurate evaluation of the two-pion contribution to 
the muon anomaly, taking into account the correlations 
between different points, is available from the refs. [13, 
|4]| . As a result, this choice guarantees a very precise 
input. We must emphasize that, once the input (H))-© 
is adopted, the treatment is optimal and no information 
is lost. A posteriori, it turns out that the results given 
by this choice are quite stringent. 

In addition to the above input from the timelike axis, 
we include the values of F{t) measured experimentally 
at some spacelike points 



P(t„) = /„ ± Sfn 



tn < 0, 



1,....,N, (8) 



where we use the most recent high precision experimen- 
tal information from [2^[23|. Thus, we will be employing 
as input Eqs. ([I])-® in order to obtain correlations be- 
tween the coefficients of the Taylor expansion ([3]). We 
will investigate also the issue of the possible zeros of the 
form factor, deriving regions where zeros are forbidden. 

In Sec. |TT]we briefly review the mathematical method 
and in Sec. IIIII the experimental information that goes 
into our computation. In Sec. IIVI we present our results 
for the parameters (c, d) and compare them with results 
available in the literature. In Sec. |V]we derive regions 
where zeros are excluded along the real axis and in the 
complex t-plane, and in Sec. I VII some discussions and 
our conclusions are presented. 



II. BASIC FORMULAE 

For solving the problem we follow a mathematical 
method presented in [6ll, . We first define the Omnes 



3 



function 



0{t) 



exp - 



dt' 



m 
t'{t' - t) 



(9) 



where 5{t) — 5\{t) for t < ijn, and is an arbitrary func- 
tion, sufficiently smooth [i.e. Lipschitz continuous) for 
t > tin. As shown in [68|], the results do not depend 
on the choice of the function S{t) for t > tin- A crucial 
remark is that the function h(t) defined by 



F{t) = 0{t)h{t) 



(10) 



is analytic in the i-plane cut only for t > The equality 
([5]), written in terms of h{t) as 



1 



dtp{t)\0{t)\^\hit)\'^a 



tin 



-TTTT 



(11) 



can be expressed in a canonical form, if we perform the 
conformal transformation 



m = 



^in \/^in ^ 
tin 's/i'm ^ t 



(12) 



which maps the complex t-plane cut for t > tin onto the 
unit disk \z\ < 1 in the z-plane defined by z = z{t), and 
define a function g{z) analytic in \z\ < 1 by 



giz)^wiz)Lu{z)F{iiz))[Oit{z))Y 



(13) 



In this relation t(z) is the inverse of z = z{t), for z{t) as 
defined in IT^ , and the last two factors give the function 
h{t{z)) defined in (fTU|) . which is analytic in \z\ < 1. Fi- 
nally, w{z) and ijj{z) are outer functions, i.e. functions 
analytic and without zeros in \z\ < 1, defined in terms 
of their modulus on the boundary, related to a/ p{t) and 
10(^)1, respectively. Equivalent integral representations 
of the outer functions in terms of their modulus can be 
written either in the z or t variables. In particular, we 



w{z) — exp 
where 

and 

a;(z) = exp 



\n\w{Q)\ 



C = exp(i6 



k(C)P = p(t(C)) 



d?(C) 



\J tin - i{z) 



\n\0[t')\&t' 



(14) 



(15) 



t.„ ^/V^^{t' -t{z)) 



Then (fTTj) can be written as 



27r Jo 



2 _ ^TTTT 



(16) 
(17) 



^From (|12p it follows that the origin i = of the t-plane is 
mapped onto the origin z = of the z-plane. Therefore, 
from (fT3|) it follows that each coefficient ^fe e i? of the 
expansion 



9{z) =90+ giz + 322^ + gj.z'^ + 



(18) 



is expressed in terms of the coefficients of order lower or 
equal to fc, of the Taylor expansion Q. Moreover, the 
values F{tn) of the form factor at a set of real points 
t„ < 0, 71 = 1, 2, N, lead to the values 

g{Zn) = w{Zn)uj{Zn) F{tn) [0(i„)]"\ Z„ = z(<„). 

(19) 

Then the norm condition (|17p implies the determi- 
nantal inequality (for a proof and older references see 

My- 



I 

6 



z 



1 



Z 



{ZIZ2)^ 
1 - ZlZ2 



6 

1 - Z1Z2 

l-z2 



(ziZat)-^ (z2ZAr) 



K 



(ziZn)'' 
1 - ZlZN 
{Z2ZN)'' 
1 - Z2ZN 



^2K 



1 - ziZn 1 - Z2ZN 



'AT 



>0, (20) 



where if > 1 is an arbitrary integer and 



fc=0 



K-l 



= g{Zn) - 9k^n- (21) 



k=0 



The same relation ([20]) holds if we replace a^'^ by an 
upper bound of this quantity and the equality sign in (jl7p 
by the < sign. Moreover, as shown in [6^, the results 
depend in a monotonic way on the value of the rhs of 
P7p . becoming weaker when this value is increased. 

The extension to the case of complex points i„, which 
enters in pairs since F(t*) = F*{t), is straightforward 
and will be discussed in Sec. |Vl 



III. EXPERIMENTAL INPUT 



We take -y/ti^ = 0.917 GeV, which corresponds to the 
first important inelastic threshold, due to the wtt pair. 
The choice of a lower value of tin is legitimate in the 
present formalism, and we will work also with ^/ti^ = 
0.8 GeV, which will allow us to compare the constraining 
power of the input conditions (|3]) and (O . 

Very precise parametrizations of the phase-shift 6l are 
given in [H, [26 1. We use as phenomenological input the 
phase parametrized as [1^ 



cotSl{t) = y^{M'^^ 



Bq + Bi 



Vt- 



t 



(22) 
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contribution to the muon anomaly for en- TABLE II: Spacelilce data from [23, [23 

ergies above 





a7 


0.800 GeV 


95.23 X 10"^" 


0.917 GeV 


22.17 X 10"^" 



t 


Value [GeV^] 


F{t) 


ti 

t2 


-1.60 
-2.45 


0.243 ± 0.012l^ go^ 
0.167 ±0.010t|; gJ? 



where k^^ = ^t/4 - and 

y/Ui = 1.05 GeV, Mp = 773.6 ± 0.9 MeV, 
Bo = 1.055 ±0.011, Bi = 0.15 ±0.05. (23) 

The function 6l obtained from is practically identical 
with the phase-shift obtained in [24;] from Roy equations 
for Vi < 0.8 GeV. The uncertainty is very small and we 
have checked that the results are practically insensitive 
to the variation of the phase-shift within the errors. 

Above tin we use in ^ a smooth phase S{t), which ap- 
proaches asymptotically tt. As shown in 68], the depen- 
dence on S{t) of the functions O and uj, defined in ((SJ and 
(IT5)l respectively, exactly compensate each other, leading 
to results fully independent of the unknown phase in the 
inelastic region. 

The two-pion contribution to muon anomaly was eval- 
uated to great precision in [40, 4lJ . The most recent eval- 
uation [4l| , based on all the available experimental data, 
gives for the total tt+tt" contribution to muon anomaly 
the value a^'' = (507.80 ± 1.22 ± 2.50 ± 0.56) x lO"!". 
In our method we need the specific contribution aJJ' of 
the energies from -y/thi to infinity. The values given be- 
low""^ are based on the B AB AR data [sj , whose spectrum 
extends up to 3 GeV. 

For the interval 0.917 - 3 GeV the two-pion contribu- 
tion is (21.73±0.24) x 10^^". Increasing the central value 
by the error, and adding an estimate of about 0.2 x 10~^° 
for the interval from 3 GeV to oo, gives the close upper 
bound a^'^ < 22.17 x 10"^° for the two -pion contribu- 
tion from 0.917 GeV to oo. As mentioned above, if we 
use in pT|) . instead of the exact value of aj^^ an upper 
bound on this quantity, the results are still valid, but are 
weaker. In order to obtain results which are in the same 
time unbiased and stringent, we need a conservative and 
accurate estimate of d'",'^ . 

For the interval 0.8-3 GeV the two-pion contribution 
in [m is (94.25 ±0.77) x 10"^". Increasing as before the 
central value by the error, and adding 0.2 x 10"^° for the 
interval from 3 GeV to oo, we obtain a^'^ < 95.23 x 10"^° 
for the two-pio n contribution from 0.8 GeV to oo. The 
final numbers for the two choices of iin are compiled in 
Table [Tj 



^ We are grateful to Bogdan Malaescu for providing us these num- 
bers. 



Finally, we use additional spacelike data coming from 
[2^ [23} . which are collected for completeness in Table 
im where the first error is statistical and the second is 
systematical. 



IV. ALLOWED DOMAIN IN THE c-d PLANE 

In this section, we present the constraints on the coef- 
ficients c and d entering the Taylor expansion ([3]) using 
the formalism developed in Sec. [ll| We list in Table 
mil the various quantities required in the basic inequality 
(|^n| , for two choices of tin ■ We implemented the normal- 
ization F(0) = 1, but kept arbitrary the charge radius 
(rj) and the spacelike values Fi and F2. Using the in- 
put from Tables [J and lllli one obtains easily from ([20)1 
a convex quadratic condition for the coefficients c and d, 
represented as the interior of an ellipse in the c~d plane. 

We consider first the constraints obtained without any 
information at spacelike points, when the determinant 
([^ has only one element, /, and the condition (|^ be- 
comes 

5o+5?+5l + 53 + --- <«r- (24) 

The quantities gi are calculated for tin — (0.8 GeV)^ us- 
ing the first line of Table [J and the first column of Table 
mil and for tin = (0.917 GeV)^ using the quantities writ- 
ten in the second line of Table [J and the second column 
of Table [mj 

In order to investigate the influence of the choice of 
the threshold tin, we show in Fig|T]the domains obtained 
with the two values of tin considered in Tables [I| and 
nil For convenience, we take (rj) = 0.43 fm^ [251 [ssl- 
53- The figure shows that the ellipse corresponding to 
tin — (0.917 GeV)^ is smaller and lies fully inside that of 
the ellipse with tin = (0.8 GeV)^, proving that the best 
results are obtained by exploiting the known phase along 
the whole elasticity region. Therefore, in what follows we 
shall adopt the choice tin = (0.917 GeV)^. 

A precise estimate (rj) = (0.435 ± 0.005) fni^ is given 
in 25]. In Fig [51 we present the domains described by 
(|20| for tin = (0.917 GeV)^ and two values of the charge 
radius (rj) — 0.43 fm^ and (rj) = 0.44 fm^ resulting from 
this estimate. The allowed domain is quite sensitive to 
the variation of (rj), being shifted towards the upper 
right end if (rj) is increased. To account for the uncer- 
tainty of the charge radius, we take as allowed domain 



5 



TABLE III: Tabulation of the quantities entering as input in H20p for obtaining the constraints on the c, d coefficients, for two 
choices of tin- The numbers Zn = z{t„) are obtained using (jf 2|) and tn given in Table Ull The numerical coefficients include the 
information on the phase below tin and the normalization F{0) — 1, while the charge radius (r^) (expressed in fm^) and the 
values Fn = F{tn) are left arbitrary. 



Quantity 


tin ~ (0.8 GeV) 


tin ~ (0.917 GeV) 


50 


0.2284 X 10"" 


0.1238 X 10""* 


51 


(0.2503(rJ) -0.0414) x 10"^ 


(0.1783(rJ) - 0.0431) x 10"^ 


52 


(0.1497c - 0.9547{rJ) - 0.II60) x 10"^ 


(0.1401c - 0.9773(r2) - 0.0985) x 10"^ 


53 


(-0.8704c + 0.3833d + 0.3879(r5) - 0.7260) x 10"^ 


(-1.0481c + 0.4712d + 0.3589{r2 ) - 0.9154) x 10"^ 


Zl 


-0.3033 


-0.2603 


Z2 


-0.3745 


-0.3285 


9{zi) 


Fi X 0.3051 X 10"'* 


Fi X 0.2066 X 10~" 


9(^2) 


F2 X 0.3984 X 10"" 


F2 X 0.2210 X 10"" 



the union of the two ellipses in Figl21 which leads to the 
ranges 

3.48 GeV"'' < c < 3.98 GeV"'', 
9.36 GeV"*^ <d< 10.46 GeV-^ (25) 

with a strong correlation between the values of c and d. 




I r:^ \ > \ I \ > \ > I 

3 3.2 3.4 3.6 3.8 4 

c [GeV"] 



FIG. 1: Comparison of the c — d domain obtained with tin ~ 
(0.917 GeV)2 and ti„ = (0.8 GeV)^ for (rj) = 0.43 fm^ 

We implement now the value at a point on the space- 
like axis, using the input given in Table [III In this case 
the determinant in (j20p has two rows and two columns. 
We choose the input at the spacelike point ti given in 
Table |TT] and account for the errors by varying Fi inside 
the error bars. In Fig. [3] we present the allowed do- 
main in the c — d plane obtained for (r^) — 0.43 fm^ and 
three values of Fi: the central value 0.243 given in Ta- 
ble [Til and the extreme values 0.265 (0.228) obtained by 
adding (subtracting) the corresponding errors added in 
quadrature. The additional information on the spacelike 
axis improves in a dramatic way the constraints on the 
c and d coefficients. The small ellipses are entirely in- 
cluded in the larger ellipse obtained without information 



10.5 — <r^> . 0.43 Im'^ 
— <r2>- 0.44 fm^ 




3.5 3.6 3.7 3.8 3.9 

c [GeV*] 



FIG. 2: Allowed domain in the c — d plane obtained with 
tin = (0.917 GeV)^ for {rl) = 0.43 fm^ and (rl) = 0.44 fm^ 

on the spacelike axis, which confirms the consistency of 
the various pieces of the input information. Varying Fi 
inside the error bars, we obtain the allowed domain of 
the c and d parameters at the present level of knowledge 
as the union of the three small ellipses in Fig. [3] This 
gives, for (rl) = (0.435 ± 0.005) fm^ the allowed ranges 

3.75 GeV"-* < c < 3.98 GeV"^, 
9.91 GeV"^ <d< 10.45 GeV"^ (26) 

with a strong correlation between the two coefficients. 
The comparison with shows that the information at 
the spacelike point improves the lower bounds on both c 
and d, a feature seen actually from Fig. [31 

Similar results are obtained using as input the second 
Huber datum t2 in Table [III Note that the formalism al- 
lows the simultaneous inclusion of several spacelike points 
in the determinant ((^0)) . In practice, as discussed in 
[13, ll^l , when more points are included the results are 
extremely sensitive to the values used as input, which re- 
quires adequate numerical methods for treating the prob- 
lem. 
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— no spacelike datum 

— F(t,) (cent) 

— F(t,) (max) 

— F(t,)(min) 




c [GeV 



FIG. 3: Allowed domain in the c — d plane calculated with 
tin = (0.917 GeV)^ and (rl) = 0.43 fm^ for three values of 
F{ti) at the spacelike point ti = —1.6 GeV^ (central value 
in Table [ll] and the extreme values obtained from the error 
intervals). Also shown is the large ellipse when no spacelike 
datum is included. 



To illustrate the issues that arise in this context, a 
small digression is in order: the formalism presented in 
this work can be used to obtain limits on the value -^(^2) 
at the second spacelike point, given the value F{ti) at the 
first one. This range results from the general inequality 
(PO)) . written as 



-9o-9i 
6 



6 



1 - 21^2 



6 

{z,Z2f 
1 - 21^2 

l-z| 



> 0, (27) 



where Zi = z{ti) and = g{zi) - go- 9iZi- 

Using as input the coefficients given in the second col- 
umn of Table HID we obtain from the above inequality 
a strong correlation between the values F{ti) and F{t2). 
For instance, taking the radius to be 0.435 fm^ and F{ti) 
at its central value in Table HIl ([27]) restricts F{t2) to the 
narrow range (0.159,0.173). The central experimental 
value of F{t2) quoted in Table |TT] is contained in this 
range, which means that the central Huber values are 
consistent with each other in the analytic framework that 
we have adopted. On the other hand, taking F{ti) at the 
lower end of the experimental interval, we obtain the al- 
lowed range of ^"(^2) as (0.135,0.153), below the experi- 
mental interval, while fixing F{ti) at the upper end yields 
the range (0.199, 0.201), above the experimental interval. 
It follows that an allowed range of F{t2) consistent with 
the experiment can be obtained only by reducing the in- 
put range of F{ti). By varying simultaneously the value 
of F{ti) and the radius, {rl) = (0.435 ± 0.005) fm^ we 
obtain for F{t2) the range (0.130,0.201), which may be 
expressed as F{t2) = 0.166 ± 0.036. The result is consis- 
tent with the numbers in Table |lll but the error is a bit 
larger than the actual experimental one. 



The digression above shows also that, by imposing si- 
multaneously the experimental values at ti and t2, we 
can only obtain a slight improvement of the allowed do- 
main of the parameters c and d. The reason is the fact 
that, as already noted above, the information on F{t2) 
forces F{ti) to lie within a slightly smaller range around 
the central value. Since the gain is expected to be small, 
we keep for simplicity as input only one spacelike point, 
which is sufficient to produce the narrow ranges reported 
in (EH). 

It is of interest to compare our predictions with pre- 
vious determinations available in the literature. First, 
the range of c given in considerably iniproves the 
bounds obtained with similar techniques in [6ll - l63l . l66j . 
The improvement is due mainly to the very accurate in- 
formation available now on the modulus, expressed in 
the values in Table ID On the theoretical side, a fit 
based on ChPT to two-loop accuracy for r decays gives 
c = (3.2 ± 0.5e:rp ± 0.9t/.eor) GcV"'' jj] . Subsequent 
calculations of the electromagnetic form factors in two- 
loop ChPT lead to the values c = (3.85 ± 0.60) GeV"" 
0, in agreement with the range given in (j26|) . and 
c = (4.49 ± 0.28) GeV"'' Q, slightly above that range. 
Finally, both the prediction c = (4.00 ± 0.50) GeV"'', 
based on the quark-mass dependence of the form fac- 
tor [6^, and the range c = (4 ± 2) GeV""* quoted as a 
conservative next-to-next-to-leading ChPT result in the 
same reference, are consistent with (|26p . On the other 
hand, a recent lattice calculation with chiral extrapola- 
tion based on two-loop ChPT gives a slightly lower value, 
c = 3.22(17) (36) GcV""* [7]. It must be noted however, 
that the lattice data are generated at rather high space- 
like momenta, t £ (—0.3,-1.7) GeV^. Therefore the 
extraction of the radius and the curvature can not be 
very precise and the corresponding uncertainties might 
be larger than estimated. 

Other determinations of the curvature are based on fits 
of experimental data with specific analytic parametriza- 
tions of the form factor. The value c = (3.90 ± 
0.10) GeV"'' was obtained in [H^ by a usual disper- 
sion relation, while a fit of the ALEPH data [l^l on 
the hadronic r decay rate with a Gounaris-Sakurai for- 
mula for the form factor [58| gives c = (3.2±1.0) GeV""*. 
Several analyses are based on phase (Omnes-type) rep- 
resentations, with various parametrizations of the phase 
along the whole unitarity cut. Their predictions, like 
c = (3.79 ± 0.04) GeV"" 0, c = (3.84 ± 0.02) GeV"" 
[ii] and more recently c = (3.75 ± 0.33) GeV"" are 
in overall agreement with ((26|) . We note also that the 
value c — (3.30 ± 0.03stat ± 0.33syst) GeV"", obtained 
recentl y fr om a fit of spacelike data with Pade approxi- 
mants [6J1 , is below our prediction f^M . It may be worth 
investigating whether the fact that the unitarity cut and 
the precise data available along it are not included in this 
analysis is responsible for the mismatch. 

As in the case of c, the range of d given in consid- 
erably improves the bounds obtained with similar tech- 
niques in (6ll - [63l [66j . The information available in the 
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literature on the cubic term in the Taylor expansion is 
not rich. Theoretical results from ChPT and lattice cal- 
culations are not yet available. From fits of the data, the 
value d = (9.70 ±0.40) GeV~^ was obtained by means of 
usual dispersion relations in [59i\ , while the Taylor expan- 
sion of the Gounaris-Sakurai parametrization (stI ]. men- 
tioned above, leads to d = 9.80 GeV~^. Both values are 
consistent with the range 



V. DOMAIN WHERE ZEROS ARE EXCLUDED 

As we discussed in the Introduction (see also [70|'). 
the formalism developed in Sec|lT] allows one to find rig- 
orously the domain where the form factor cannot have 
zeros. The method amounts to testing the consistency 
of the assumption that a zero is present with the other 
pieces of the input. Let us assume first that F(t) van- 
ishes at some real point to. From ([T^ it follows that 
g{zo) — 0, where zq — z{ta). We therefore include this 
value in the determinant (pil)) and test the validity of the 
inequality: if it is satisfied, a zero is possible, if it is vio- 
lated, the zero is forbidden. In particular, if we use only 
the information on the normalization F{0) = 1 and the 
charge radius, with no input on the spacelike axis, we 
obtain from ([20)1 and (|2ip the following condition 



-go -.912=0 



-go-gizo 



1 



<o. 



(28) 



for the points zq such that the form factor cannot vanish 
at to = i{zo). Here go and gi are expressed cf. Table Hill 
in terms of the charge radius. 

If we include in addition the value at a point zi = z{ti), 
the condition reads 



-go-gizo 



-go~giza 
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l-zl 


1 - ZqZi 




{zif 


1 - ZoZi 


l-zl 



< 0, (29) 



with ^1 = g{zi) - go- ffi^i- 

With the values given in Tables H] and Ull] for tin 
(0.917 GeV)2 and (r^) = 0.43 fm^ the inequality 1^ 
implies that simple zeros are excluded from the inter- 
val -1.93 GeV^ < to < 0.83 GeV^ of the real axis. If 
we impose the additional constraint at a spacelike point 
ti = —1.6 GeV^, the interval for the excluded zeros is 
much bigger. The left end of the range is actually quite 
sensitive to the input value Fi — F[ti). Using the cen- 
tral value Fi — 0.243 given in Table [TTl we find from 
that the form factor cannot have simple zeros in 
the range -5.56 GeV^ < to < 0.84 GeV^ By vary- 
ing Fi inside the error interval given in Table |ll] (with 




FIG. 4: Comparison of the domains without zeros obtained 
from pO|) using tin ~ (0.8 GeV)^ (smaller domain) and tin = 
(0.917 GeV)^ (bigger domain), for {rl) = 0.43 fm^. 



errors added in quadrature), we find that zeros are ex- 
cluded from the range -12.67 GeV^ < to < 0.84 GeV^ 
for Fi = 0.265 at the upper limit of the error inter- 
val, while for the lower limit Fi — 0.228 the range is 
-4.46 GeV^ < to < 0.84 GeV^. 



E 




FIG. 5: Domain without zeros obtained from H30|) using 
tin = (0.917 GeV)^, for two values of the pion charge radius, 
(rl) = 0.43 fm^ (smaller domain) and (rl) = 0.44 fm^ (bigger 
domain) . 



We now turn to the study of complex zeros. The for- 
malism presented in Sec. |lT]can be easily adapted to in- 
clude complex values of the form factor outside the real 
axis. Since the form factor is real analytic, its zeros occur 
in complex conjugate pairs, i.e. if F(to) = 0, then also 
F{to) = (a double zero occurs as to approaches the real 
axis) . One can show that the determinant condition 
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E 




FIG. 6: Domain 
(0.917 GeV)^ and {rl) = 0.43 fm^ 
tral experimental value F(ti) 
ti = -1.6 GeV^ 



without zeros obtained with tin ~ 
using in addition the cen- 
0.243 at the spacelike point 




-20 -10 
Ret 



FIG. 7: Comparison of the domains with no zeros obtained 
with tin = (0.917 GeV)^ and (rl) = 0.43 fm^, for the spacelike 
input Fi = 0.265 (bigger domain) and Fi = 0.228 (smaller 
domain) . 



for the domain without zeros is generalized to 

^2 „2 



~9o~9iZo 
-go-gizo 



-.90 -.91^0 
-A 



z, 







1-Z2 



-.90 -.91^0 

(^0*)^ 

1 - iz*oy 

1-koP 



<0. (30) 



The determinant is real since the corresponding matrix 
is Hermitian. The 4x4 determinant that includes in 
addition a value at a spacelike point ti can be easily 
written down. 

We first apply the inequality (pO| to illustrate the de- 
pendence of the domain without zeros on the value of iin 
used in the calculations. As seen from FigHl the larger 
value iin = (0.917 GeV)^ leads to a domain that extends 
to high values of \t\ in all the directions of the complex 
plane, which shows that, like in the case of the c — d do- 
main, the best results are obtained if the phase condition 
(H]) is used along the whole range of validity. 

The dependence of the domain on the variation of (r^) 
is shown in FiglS] As expected, for a larger charge radius. 



{rl) = 0.44 fm^ the zeros are excluded from a bigger 
complex domain around the timelike axis, while the left 
end of the domain, around the spacelike axis, is almost 
insensitive to the slope at i = 0. 

The effect of an additional input at a spacelike point 
is illustrated in FiglHl where we show the domain in 
the complex plane where zeros are excluded, using ti — 
— 1.6GeV^ and the value F(ti) = 0.243 (the central ex- 
perimental value given in [l^, [23|). By comparing FiglS] 
with the large domain in Fig|31 one can see that the 
knowledge of the form factor at a spacelike point excludes 
zeros in a larger domain near the spacelike axis, while it 
has a smaller influence on the right part of the domain. 
This feature is present also in Figl?! which shows the 
sensitivity of the domain to the input value of F{ti). As 
is seen in the figure, the larger value F{ti) — 0.265 ob- 
tained from the upper limit of the error bar, excludes the 
zeros in a domain extending to considerably larger values 
along the spacelike axis. 

The results on the zeros reported in the literature 
[H, m-El 113, El] are rather controversial. The best 
results for the regions free of zeros were obtained in 
[12, 113, SI], by a method related to ours. However, 
since the experimental information at that time was poor, 
the authors were forced to make some ad-hoc assump- 
tions, especially on the modulus on the timelike axis. At 
present the precise measurement of the modulus gives a 
solid basis to our results. 

The issue of zeros is of relevance for the analytic rep- 
resentation of the form factor using phase (Omnes)- or 
modulus- type representations, which require the knowl- 
edge of the zeros. Such representations were extensively 
studied in the past JJ, 46, 50, 52, 53, 65j], and often 
are based on the assumption that zeros are absent. Our 
results, which show that the zeros are excluded from 
a rather large region at low energies, give support to 
such representations, and confirm also theoretical expec- 
tations based on ChPT or more general physical argu- 
ments [Hi. 



VI. DISCUSSIONS AND CONCLUSION 

The experimental information available at present 
on the pion electromagnetic form factor is very rich. 
The recent high statistics measurements of the modulus 
by BABAR and KLOE collaborations [s^-lli], supple- 
mented by the phase in the elastic region known with 
accuracy from the P-wave of tttt scattering [23 - [26l |. are 
expected to considerably constrain the behavior on the 
timelike axis. The values of the form factor on the space- 
like axis are also measured with increasing precison. The- 
oretically, predictions on the pion form factor at low en- 
ergies are available from ChPT and lattice QCD, while 
perturbative QCD predicts the behavior at high energies 
along the spacelike axis. The transition to the perturba- 
tive regime is known to be an open problem that deserves 
further study in the case of the pion form factor. 
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Analyticity is the ideal tool for connecting the low- and 
high-energy regimes for physical quantities like the pion 
form factor. The full treatment of the present rich experi- 
mental and theoretical input, which might overconstrain 
the system, is a challenge for the future investigations 
based on analyticity. In the present work we do not per- 
form such a complete analysis, but exploit only in part 
the present information on the modulus on the unitarity 
cut. However, even in this limited frame we obtain quite 
stringent conclusions on the low energy properties of the 
form factor. 

The conditions used as input in our approach are ex- 
pressed by the phase condition ^ and the integral of 
the modulus squared ([51), which we further restricted by 
choosing the weight p{t) as the kernel relevant for the 
two-pion contribution to the muon anomaly, cf. ^ and 
([7]) . A more general class of suitable weights will be in- 
vestigated in a future work. Once the input is chosen, it 
is exploited in an optimal way by a mathematical formal- 
ism presented in Sec. |TT1 leading to strong correlations 
between the coefficients of the Taylor expansion (jSj at 
t = and the values of the form factor on the spacelike 
axis. 

Our basic results are contained in Eqs. (|20| and (|2T]) . 
where the input quantities are defined in Tables Hl lIIII 
The numerical coefficients in Table IIIII depend on the 
normalization -F(O) = 1 and the phase of the form factor 
below the inelastic threshold tin, being vary stable with 
respect to small variations of the phase. Moreover, as 
emphasized in Sec. [TTl the results are independent of 
the unknown phase of the form factor above the inelastic 
threshold tin- In Table Hill the charge radius (r^), the 
higher-order Taylor coefficients c and d, and the values 
of the form factor at several spacelike points are kept free, 
so the formalism can be easily applied for finding model 
independent correlations between the values of the form 
factor at different points and for testing the consistency 
of input values known from different sources. 

In Sec. IIVI we derived stringent constraints on the 
allowed values of the higher-order coefficients c and d 
of the Taylor expansion The best results are ob- 

tained with tin = (0.917 GeV)^, which corresponds to 
the physical inelastic threshold produced by the lot: chan- 
nel. The charge radius and an additional information 
at a spacelike point were used as input. In (|25p and 
([26l) and in FigslT]- Fig. [3] we illustrated the results for 
{rl) in the range (0.43 - 0.44) fm^ and F{-1.6 GeV^) = 
0.243 ± 0.012to°ll [H Hi- It is remarkable that the 
allowed ranges are already comparable in precision with 
other determinations in the literature based on specific 
parametrizations . 

The present method can be used also to obtain bounds 
on the values of the form factor along the spacelike axis, 
using as input the information on the timelike axis, to- 
gether with some values inside the analyticity domain. 
As discussed in Sec. |IVl using as input the value F{ti) 
at the first Huber point, we obtain stringent limits on the 
value F{t2) at the second point, with a strong correlation 



between the two. Of course, the method can be applied 
in principle also to higher spacelike energies. However, 
with our choice of the weight p, we expect that the pre- 
dictions will become gradually weaker when the energy is 
increased. Indeed, since p decreases rapidly at large mo- 
menta, the condition ([5]) provides stringent constraints 
on the low energy parameters like c and d, but in the 
same time it imposes weak constraints on the behavior 
of the form factor at large energies. A different choice 
of p could lead to interesting results also for the behav- 
ior at higher energies, but this is beyond the scope of 
the present analysis and will be investigated in a future 
work. 

In Sec. |V] we showed that the same formalism leads 
to an analytic description for the regions of the complex 
plane where the zeros of the form factor are forbidden. 
Our results are contained in Eqs. (l^5)) - (l5ni) and are il- 
lustrated in Figs. dH) - d?]), for the same input (rj) and 
F{ti). We obtain a rather large domain where zeros are 
excluded, which gives support to Omnes-type representa- 
tions, which often assume the absence of the zeros. Our 
results also confirm theoretical expectations on the ab- 
sence of zeros at low energy, based on ChPT or general 
physical arguments [11] . We note that by our method we 
can find rigorously the domains free of zeros, but we can 
say nothing about the remaining domains, where zeros 
may be present or absent. Alternative methods, based 
on modulus representations [H, Ull , can rule out in 
principle the zeros from the whole complex plane pro- 
vided they are absent. However, these methods are very 
sensitive to the input and led to controversial results in 
the past. An update of such analyses using the recent 
precise determination of the modulus would be of much 
interest. 

We finally note that the mathematical formalism ap- 
plied in this paper may be useful also for finding an an- 
alytic parametrization of the form factor suitable for fit- 
ting the rich amount of experimental data. Namely, the 
representation of F{t) that results from ([T^ involves the 
known functions ^(z), which accounts for the weight p{t), 
Lu{z) and 0{t), which implement the phase below tin, and 
the arbitrary function g{z), analytic in the t-plane cut 
for t > tin, or equivalently in the unit disc \z\ < 1 of 
the z-plane defined by the conformal mapping (|12p . The 
expansion (|18p is convergent in \z\ < 1, and moreover the 
coefficients satisfy the inequality ([Ml) , which is very use- 
ful in order to control the higher orders of the expansion 
and the truncation error. 
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